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Testing a Topology Conserving Gauge Action in Lattice QCD 
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We study lattice QCD with a gauge action, which suppresses small plaquette values. Thus the MC history 
is confined to a single topological sector over a significant time, while other observables are decorrelated. This 
enables the cumulation of statistics with a specific topological charge, which is needed for simulations of QCD in 
the e-regime. The same action may also be useful for simulations with dynamical quarks. The update is performed 
with a local HMC algorithm. 



1. MOTIVATION 

Our goal is to explore the applicability of lat- 
tice gauge actions for QCD, which are designed 
such that small plaquette values are strongly sup- 
pressed. If such an action can be identified, we 
expect the following virtues: 

• an acceleration of dynamical fermion simu- 
lations 

• control over the topological charge . 

In particular we hope for applications in the 
e-regime of QCD p^. In that regime, the pion 
Compton wave length clearly exceeds the box 
length, m^ 1 3> L, which is an unphysical situ- 
ation. However, simulations in a such small vol- 
ume may provide physically significant informa- 
tion, since the low energy constants of chiral per- 
turbation theory take the same values as in a large 
volume. Hence the hope is to evaluate them with- 
out the requirement of a large lattice. 

However, this implies several conditions for the 
lattice formulation: the lattice fermions should 
keep track of the chiral symmetry and give ac- 
cess to very light pions. Moreover, the topologi- 
cal charge Qtop needs a sound definition, since the 
measurements are performed at distinct Qtopl [3 

Both of these requirements are provided by 
Ginsparg- Wilson fermions: they have an exact, 
lattice modified chiral symmetry 3 k , and the 
fermionic index defines Qtop 0]- Their simula- 
tion is now possible, at least quenched [1] . Recent 
results were obtained for the Dirac spectrum [Hj 
and for meson correlation functions 0. 



In particular the Neuberger Dirac operator [H] 
can be applied. Its index — and hence the topo- 
logical sector — cannot change as long as each 
plaquette variable Up obeys the constraint |9I10| 



S P := 1 - -Re Ti(U P ) < e 



1/20.49 



where (3Sp is the plaquette action. 

Therefore Qtop is fixed under continuous defor- 
mations of the gauge configuration if we use the 
plaquette action /3S a , with 



S a (Up) 



Sp(Uf 



[1-Sp {U P )/e\ 
+0O 



< e 



S P (U f 



otherwise 



(1) 



for a > 0. This lattice gauge action (for a = 
1) was introduced by M. Liischer for conceptual 
purposes and applied by Fukaya and Onogi 
in Schwinger model simulations |13j . 

The use of such a lattice gauge action has the 
advantages that the continuum property of stable 
topologies is reproduced, and that tedious com- 
putations of the index can be saved. Finally it 
allows for the cumulation of statistics in a spe- 
cific topological sector; first experience in the e- 
regime shows that in particular | Qtop 1 = 1,2 are 
useful (the sector Qtop = suffers from strong 
fluctuations, and at |Qtop| > 2 quenched chiral 
perturbation theory |12| fails in volumes with box 
length L< 1.5 fm). 

On the other hand, an obvious problem is that 
due to the constraint implemented in S a (Up) the 
physical lattice spacing tends to be very small. 
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Figure 1. Above: ratio between the standard pla- 
quette action Sp and the modified plaquette ac- 
tion S a . Below: ratio between forces in configura- 
tion space, which drive the local HMC algorithm. 

For practical simulations we hope that e values 
clearly above the theoretical bound still suppress 
topological transitions sufficiently. For instance, 
in the Schwinger model e = 1 already stabilized 
Qtop over hundreds of configurations [T3*| . 

2. A LOCAL HYBRID MONTE CARLO 
ALGORITHM 

Since gauge actions of the type Q are non- 
linear in the link variables, the heat bath algo- 
rithm cannot be applied. Instead we use a local 
HMC algorithm Compared to the Wilson 

action, the force is just changed by an extra fac- 
tor on each plaquette, 

= SS a (U P ) = 6Sp(U P ) I + ^Sp 
SU SllA SU Xtfl ' (1 - Sp/e)^ ■ 

The new plaquette action S a and the force F a are 
illustrated as functions of the standard plaquette 
action Sp in Fig. Q 

3. RESULTS 

We set a = 1 and as a first experiment we 
searched for the line of a constant plaquette vari- 
able (Sp) on a 4 4 lattice, as (3 and the action 
parameter e are varied. The result is shown in 
Fig. |21 As we decrease e, very small values of (3 
are needed to keep the plaquette constant. 
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Figure 2. The lines of constant plaquette values 
in the plane spanned by 1/e and (3, on a 4 4 lattice. 
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Table 1 

Results for the physical scale ro/a on a 16 4 lattice 
at different values of 1/e and (3. We add the cor- 
responding [3 value for the Wilson gauge action, 
which we denote as (3w ■ Finally we display the 
autocorrelation time with respect to IQtop 1 !; an< ^ 
with respect to the plaquettes, in our MC history. 

As a more serious approach to identify a line of 
a constant physical scale, we now proceed to a 16 4 
lattice and measure r /a (at r = 0.5 fm) follow- 
ing the standard procedure, see e.g. Ref. [El- Fi- 
nite size effects are expected to be on the percent 
level for our results presented in Table (though 
we only indicate the statistical error). 

We also performed first tests of the topological 
stability. In order to arrive at a quick first im- 
pression this was done with cooling and searching 
for the first plateau of the action (which is not 
sensitive to the sign of Qtop)- The histories for 
IQtop 1 ! — ignoring instanton/anti-instanton can- 
cellations — are shown in Fig. |3| 

We recognize a clear progress in view of the 
topological stability as s decreases along a line of 
approximately constant physics. This can also be 
quantified by the autocorrelation in these histo- 
ries, as the fifth column in Table shows. 

On the other hand, the last column in Table 
H shows that the decorrelation with respect to 
the plaquette values becomes even better as e de- 
creases (at a fixed scale). 
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Figure 3. Histories of \Q c t °° l \ on a 16 4 lattice 
for the sets of parameters that appear in Table 
1. They correspond approximately to the same 
physical scale. 

4. CONCLUSIONS 

"Topology conserving gauge actions" could be 
highly profitable in QCD simulations. The sup- 
pression of small plaquette values may speed up 
the simulations with dynamical quarks. A stable 
Qtop is useful in particular in the e-regime. 

We are exploring the applicability of such ac- 
tions, in view of the physical scale and the topo- 
logical stability. This is an ongoing project. 
At present we have first promising candidates 
for suitable parameters which provide topological 
stability to some extent, while keeping the phys- 
ical lattice spacing at a reasonable value. 

Future tests will involve index measurements to 
verify the topological stability. At last we men- 
tion that the gauge action JJ) has the problem 
that once a plaquette violates the constraint, a 
force in the wrong direction sets in. Hence such 
configurations had to be rejected. This happens 



more frequently as e is decreased. We now want 
to test variants of the action QJ, where the de- 
nominator is replaced by a quadratic or exponen- 
tial factor in order to avoid this problem, and also 
in view of the caveat pointed out in Ref JH] . 

We thank Martin Liischer for helpful comments, 
and the Deutsche Forschungsgemeinschaft for finan- 
cial support through SFB-TR 9-03. 
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